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' ^-statistics is defined in the context of the Lie algebra sl(n + 1). Some thermal 

\0 ■ t. V / 

l> . 

^SJ , teract only via statistical interaction imposed by the Pauli principle of j4-statistics. 

' Apart from the general case, three particular examples are studied in more detail : (a) 

the particles have one and the same energy and chemical potential; (b) equidistant en- 
. ergy spectrum; (c) two species of particles with one and the same energy and chemical 

I potential within each class. The grand partition functions and the average number of 

particles are among the thermodynamical quantities written down explicitly. 
Oh! 



1 Introduction 

The first attempts to generalize canonical quantum statistics go back to Gentile |l]], who 
considered statistics, which are intermediate between Fermi-Dirac (FD) and Bose-Einstein 
(BE) statistics. More precisely. Gentile introduced statistics with the property that the 
maximal occupation number of particles on any orbital is larger than 1 (hence the statistics 
is not FD), but is finite (hence the statistics is not BE). Since that time various generaliza- 
tions of quantum statistics have been proposed both in quantum field theory 0, |^, ^ and 
in condensed matter physics ^, 0, some of them inspired by new developments in con- 
formal field theories and related lattice models ( and references therein) and in quantum 
groups 1^, |TD[. For an overview of generalized quantum statistics formulated in terms of 



deformed algebras or generalized Fock spaces, we refer to |TT], |T2 . 

In 1950, Wigner has shown (on a simple example) that there might exist statistics, 
which are compatible with the principles of quantum theory without the necessity that the 
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position and the momentum operators satisfy the canonical commutation relations. This 
more general statistics discovered by Wigner turned out to be the para-Bose statistics of one 
pair of creation and annihilation operators (CAO's) |]13|. Three years later Green introduced 
both para-Bose (pB) and para-Fermi (pF) statistics in the more general frame of quantum 
field theory 

In the present paper we study the macroscopic properties of a certain type of statistics, 
called ^-statistics. It was introduced in M, and studied further from the microscopic 



point of view in |T6|. A-statistics resembles the pF statistics insofar as the creation and 
the annihilation operators of both statistics generate simple Lie algebras : any n pairs of 
parafermions generate the orthogonal Lie algebra so{2n + 1) = i?„ fTTl , |T^, whereas any 
n pairs of A-CAO's generate the Lie algebra sl{n + 1) = (which justifies the name A- 
statistics). A-statistics resembles also Bose statistics : similar to bosons, the A-creation 
(resp. annihilation) operators commute with each other. The Fock representations for pF, 
pB and A-statistics are constructed in one and the same way : they are generated out of a 
vacuum by creation operators only. The Fock representations in all three cases are labelled 
by a positive integer p = 1,2, . . ., called the order of statistics. Moreover the metric within 
any Fock space is defined with the usual Fock space technique. It is essential to point out 
that contrary to the CAO's of parastatistics, the A-creation operators , . . ., a+ (resp. the 
A-annihilation operators a^, . . ., a~) commute with each other. For this reason (apart from 
the trivial case) they differ essentially from the CAO's of the g-ons |]19| or from the CAO's 
associated with solutions of the spectral Yang-Baxter equations [20|, since in these works 
relations of the type a~a~ = Rija~a~ are imposed. 

In the case of para-Fermi statistics of order p no more than p particles can be accommo- 
dated on any orbital. The filling of the orbitals is however completely independent of each 
other. Here comes one of the essential differences with A-statistics. The Pauli principle for 
A-statistics says that if the order of statistics is p, then the system cannot accommodate 
more than p particles. Thus, if p = 10 and 10 particles are already accommodated on the 
first orbital, then no more particles can be added to any orbital. For this reason A-statistics 
gives perhaps the simplest example of an exclusions statistics P, |^ : the number of the 
available places on a certain orbital depends on how many particles (independently where) 
are already accommodated in the system (see |1^ for more discussions of this issue). Here 



it is an appropriate place to say that the word particle is used in the context of this paper 
as a collective name for particles, quasiparticles, excitations, etc. 

In Section 2 we recall shortly the definition and the main microscopic properties of A- 
statistics. Similarly as for para-Fermi statistics the creation and the annihilation operators 
a^, . . . ,a^ of sl{n + 1) are defined via triple commutation relations (see (^?l|) ). These triple 
relations define completely the Lie algebra sl{n + 1), a property which was indicated for 
the first time by Jacobson |2l|. For this reason we call the CAO's of A-statistics Jacobson 



generators. 

In Section 3 we write down explicitly the sl{n + 1) grand partition function Z{p,n) and 
the average number of particles in the system N{p,n), see equations ( ^.9] ) and ( |3.16D , under 
the general assumption that the energy of each particle on orbital i is e^. In this context 
n is the number of orbitals of the system and p is the order of the statistics, a positive 
integer, which labels the inequivalent Fock space representations, see (p.4|). Because of the 
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Pauli principle the orbitals cannot be considered as independent subsystems (as in BE or FD 
statistics) : the filhng of any orbital depends on the states of the other orbitals. Therefore 
we derive the thermodynamical quantities directly for the n-orbital system, assuming that 
it is in a thermal and diffusive contact and in a thermal and diffusive equilibrium with a 
much bigger reservoir. As we shall see, the k-th complete symmetric functions, see ( |3.5D , 
turn out to be a particularly convenient tool for the description of the thermal properties of 
the system. 

In the remaining three sections we consider different specializations of the general set- 
tings of Section 3. First (Section 4) we assume that all orbitals (i.e. single particle states) 
have one and the same energy and chemical potential. We express the grand partition func- 
tion and the ensemble average number of particles via hypergeometric functions (see, for 
instance, equations ( |4.5| ) and ( |4.17| )). Two special cases are considered in some more detail. 
The first one corresponds to n = 1. Here, the p = 1 representation leads to the Fermi-Dirac 
distribution function, whereas p = oo corresponds to the Bose-Einstein distribution. For 
all other values of p the distribution function is intermediate between the FD and BE dis- 
tributions. The second case, corresponding to p = 1, see Figure 2, leads to the so-called 
hard-core fermions or hard-core bosons (locally they coincide). Such particles are natural 
ingredients in multi-band Hubbard or various Heisenberg spin models, where configurations 



which contain more than one particle on each lattice site are strictly prohibited (see [|T6| for 
more discussions on the topic). 

In Section 5 a model with equidistant energy levels is considered. The orbitals are labelled 
by the energy. The grand partition function is written in terms of the so-called g-generalized 
or basic hypergeometric functions, see ( |5.1CI| ). The conclusion is that for big energy gaps or 
at very low temperatures all particles "condensate" on the lowest energy orbital. The case 
with p = 1 is considered in more details. 

In Section 6 we consider two species of particles. Those of the first kind A (resp. of kind 
B) have one and the same energy (resp. e^) and chemical potential fx a (resp. /i^). Apart 
from the grand partition function and the average number of particles N{p,n), also the 
thermal average N{p, n)A of the number of particles of kind A and of kind B are computed. 
On the example of s/(5) with p = 4 the general accommodation properties are demonstrated. 
For instance the region with < fiA and > /i^ is populated most probably with particles 
of the first kind, see Figure 4, whereas the region with < f^A and €3 < fJ'B is populated 
with approximately the same number of particles of both kinds, see Figure 3. 

Throughout the paper we use the following abbreviations and notation (some of them 
standard) : 

CAO's : creation and annihilation operators; 
GPF : grand partition function; 
N : all positive integers; 
[a, b] = ab — ba. 
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2 Microscopic properties of A-statistics 



In this section we list shortly the basic definitions and some of the microscopic properties of 
A-statistics. In particular, we shall define 

• the CAO's of A-statistics and their "triple commutation" relations, 

• the Fock spaces of A-statistics, and the corresponding Pauli principle, 

• the Hamiltonian being studied in these Fock spaces. 



For more details and a derivation of the results we refer to |ll4| , |T5|, |T6|. 

The CAO's of A-statistics are equal to the Jacobson creation and annihilation operators 
a^, af, . . . , of sl{n + 1), which are defined as 2n operators satisfying the relations 

[K^ «7]' «fe ] = ^kjai + 6ija+, 

[[at, aj], ctfe ] = SkiaJ - 6ija^, (2.1) 
[a+,a+] = [a;,aj] = 0. 

The sl{n + 1) generators expressed in terms of the Jacobson CAO's read : 

eio = af, eoi = a,", e^j - Cqo = [a+, a^], aj = [a^ , aj]; i j = 1, . . . , n. (2.2) 

Above {eij\i,j = 0, 1, . . . , n} are the known Weyl generators of gl{n + 1) : 

i^ij, Cki] = SjkCii — SiiCkj- (2.3) 

As in the case of parastatistics the Fock spaces W{p, n) of A-statistics are labelled by 
an order of statistics p, where p runs over all positive integers : p G N. Each state space 
W{p,n) is defined by the requirement that it contains a vector |0), a vacuum, such that 

«r<|0) = MO), %|0) = 0; pen, i,j,k = l,...,n. (2.4) 

The Fock spaces are finite-dimensional irreducible sl{n + l)-modules. All vectors 

{at)Haty^...{aty-\0) (2.5) 

subject to the restriction 

ll + l2 + --- + ln<P (2.6) 

constitute a basis in W{p,n). 

Here a remark is in order. The linear span of all vectors ( |2.5|) for any G 
{0, 1, 2, . . .}, namely without the restriction ( p.6|) , is an infinite-dimensional s/(n+ l)-module 
W{p,n). The latter is however not irreducible. W{p,n) contains an (infinite-dimensional) 
invariant subspace Winv{p,n), which is the linear envelope of all vectors (|2.5|) with h + I2 + 
... + /„ > p. Then W{p,n) is a factor module of W{p,n) with respect to Winv{p,n) (and 
the vectors (|2.5|) subject to the restriction (|2.6|) are representatives of the corresponding 
equivalent classes in W{p,n)/Winv{p,'n)). 
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Define a Hermitian form ( , ) on W{p, n) witli tlie usual Fock space tecfinique, namely 
postulating (in addition to a~|0) = 0) that 



(a) (0|0) = 1, 

(b) (0|a+ = 0, 2 = l,...,n, 

(c) {{atr-{4r- ■ ■ ■ («+)-" |o), {attiatt " " " {<tm = 

ma-r- ■ ■ ■ {a^r^{aY)-^{at)Haty^ • • ■ (a+)'"|0)) 



, + Vn 



(2.7) 



With respect to this form any two different vectors ( p. 51 ) are orthogonal. All vectors 



\p; h,...Jn) 



■|0), /l + /2 + --- + /n <P (2. 



constitute an orthonormal basis in W{p,n), i.e. ( , ) is a scalar product. Moreover the 
Hermitian conjugate to a~ is af, {a^)* = af, which is an important physical requirement. 
The transformation of the basis (12. 81) under the action of the Jacobson CAO's reads : 



\Pt ^ly ■ ■ ■ ■> hy ■ ■ ■ 1 ^n) 



Oi \p;h,...,h,...,ln) 



\ 



(/i + l)(p-^/,) \p-h...,li_i,h + l,l ), (2.9) 



k{p + 1) b; ^1 • • • , k-i, k - 1, k+i In). (2.10) 



For further use we extend W{p,n) to an irreducible gl{n + 1) module, setting (below and 
throughout Ni = en, i = 0,1, . . . ,n) 



Nq\p;Ii,12, . . . ,ln) = {p-^li)\p;li,l2,---Jn)- 



(2.11) 



i=l 



Then 



Ni\p;li,l2,...,ln)=k\p;hj2,---Jn), i = l,...,n. (2.12) 

The basis vectors \p;li, . . . ,ln) in W{p,n) are in one to one correspondence with all 
distinct n-tuples (Zi, . . . , In) with integer non-negative entries li, . . . ,1^ such that li + - ■ ■+/„ < 
p. Based on this we often write {li, . . . , In) instead of |p; /i, . . . , 

In the present paper we will study some macroscopic properties of A-statistics for a 
Hamiltonian which is a simple sum 

n 

H = Y,^,Ni. (2.13) 

i=l 

This Hamiltonian can also be written entirely via creation and annihilation operators : 

^ n n 

H=^—lY.'^{p + n[at,aT]- J][a+a,-]). (2.14) 



i=l 



Clearly, H is an element from the Cartan subalgebra of gl{n + 1). Since H\Q) = 0, the 
energy of the vacuum is zero. The commutation relations of H with the CAO's read : 



[H,af]=±e^af. (2.15) 

If \E) is a state with energy E, then 

Haf\E) = iE±e..)af\E) (2.16) 

and therefore each af (resp. a^) can be interpreted as an operator creating (resp. annihilat- 
ing) a particle (quasiparticle, excitation) on orbital i (with energy e,). Since 

H\p; li,l2, ...,ln) = {eih + 62/2 H h eJn)\p; hj2, ■ ■ -Jn), (2.17) 

|p; /i, /2, • • • , In) is interpreted as a state with li particles on the first orbital, I2 particles on 
the second orbital and so on, /„ particles on the last orbital. 

The restriction ( |2.(j| ) expresses the Pauli principle of A-statistics in W{p, n). It says that 
the system can accommodate up to p, but no more than p particles. For this reason A- 
statistics falls into the class of exclusion statistics in the broad sense : the number of the 
allowed particles to be accommodated on a certain orbital depends on the number of the 
particles that have already been accommodated in the system. This is perhaps the simplest 
form of a statistical interaction : the Hamiltonian (12. 131) has the form of a "free" Hamiltonian 



and the interaction is introduced via a change of statistics. It will be interesting to find out 
whether one can obtain the same results adding to the Hamiltonian ( |2.13| ) an interaction 



term and changing the statistics to Bose statistics. It is known that a similar phenomenon 
can take place in quantum mechanics [53]. 



3 The grand partition function 

Here we shall study some macroscopic properties of A-statistics. For our considerations it is 
irrelevant whether the different orbitals correspond to different particles, to different energy 
levels of particles of the same kind or to different internal states of the particles. The only 
assumption is that they satisfy the Pauli principle for A-statistics. 

As usually, we assume that the system is in a thermal and diffusive contact and in a ther- 
mal and diffusive equilibrium with a much bigger reservoir. Denote by r its (fundamental) 
temperature and let /Xj be the chemical potential for the particles on orbital i. 

The general principles (and approximations) of statistical thermodynamics assert that 
the probability P(p, n; r) for the system to be in a (quantum) state r = (/i, ...,/„) with the 
number of particles Nr = h + ■ ■ ■ + In and energy Er = hei + ■ ■ ■ + is given by the 
expression : 

exp (Er=i^"H/^i - 

V{p, n; r) = ^ ^ ^. (3.1) 



The numerator in (|3.1| ) is the Gibbs factor of the system in the state r = (Zi, . . . , /„) and 



Z{p,n) is the grand partition function (GPF), namely the sum of the Gibbs factors with 
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respect to all states (/i,...,/„) of the system, i.e. over all possible non- negative integers 
li, . . . Jn such that < /i + ■■■ + /„< p : 

Z{p,n)= (exp(^^^))'Hexp(^^^))'^---(exp(^^^^))'". (3.2) 

0<h+-+l, 

In terms of the notation 



T T T 

0<h+-+ln<p 



Xi = exp ^ ^' )' i = l,---,n, (3.3) 

we rewrite ( |3.2| ) as follows, 

p 

Z{p,n)= Yl ^1^2 ■■■^n=Y. Yl ^1^2 (3.4) 

In the general setting, which we consider so far, it is appropriate to introduce the complete 
symmetric functions hk{xi, . . . , x„), k = 0,1, . . ., which play an important role in the theory 
of symmetric functions ||2^. The k-th complete symmetric function hk{xi, . . . , x„) is the sum 
of all distinct monomials of total degree k of the variables xi, X2, . . . , x„ : 

hkiXl, . . . , Xn) = ^ ^ 2^1 " ' ' . (3-5) 

For example, h^^xi, X2, x^) = 1, hi{xi,X2, X3) = xi + X2 + x^, 

h2{xi,X2, X3) = x\ + x\+ x'l + X1X2 + XiXz + X2X3. 

In terms of hk{xi, . . . , x„), the GPF Z{p, n) reads : 

p 

Z{p,n) =Y,hk{xi,. . . ,Xn). (3.6) 

fc=0 

Clearly, hk{xi, . . . ,Xn)/Z{p,n) yields the probability for the system to contain k particles. 
In order to evaluate the sum ( |3.6| ) we use the following generating function |^ (1.2.5)], 



k=0 

Now compute 



00 ^ 



00 00 p 

p=0 p=0 k=0 

00 00 00 00 

= ^ ^ hkiXi, Xn)t^ = X] hk{Xi, Xn)tP^'' 

k=0 p=k k=0 r=0 

00 00 

= (5^/ifc(a;i,...,x„)t'=)(^r) 

fc=0 r=0 



1 -j^ 00 

X2t)...{l-Xj)l-t = 5 • • • ' 
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Hence 

p 

^ ^ hki^^ly • • • ) ^n) hpi^Xi, . . . , Xfi, 1) hpi^Xi, . . . , Xi—i, 1, Xj, . . . , 3^^)- ("^•^) 
fe=0 

We have written the last term in the rhs of ( p.8|) for further use. It follows from the property 
that hp is symmetric with respect to its arguments and therefore these arguments can be 
reordered in an arbitrary way. 

Applying ( p.8| ) to ( p^ ) we obtain 

p 

Z{p, n) = ^ hk{xi, ...,Xn) = hp{xi, ...,Xn,l)- (3.9) 

Using the GPF (|3.9| ) one can determine various other thermodynamical quantities and 
in particular the average number of particles in the system. 

According to the probability V{p,n]li, . . . ,ln) for the system to be in the state 

r = [li, ...,/„) with Nr = h + ■ ■ ■ + In particles reads (in terms of the variables Xi) : 

fv»^l ... 

V{p,n;h,...,Q = ^^l^^^. (3.10) 

Then the average number of particles in the system is 

N{p,n) = J2 {h + --- + ln)V{p,n;h,...Jn) 

0<h+---+l„<p 

. . . rp^n 

= E ('■ + - + «%^. 

0<ll + - + ln<p ' 

which can also be written as 

n n 

N{p,n) = "^Xkd^^ In Z{p,n) = ^rd^^ In Z{p,n). (3.11) 



k=l k=l 



Since 



J2 {h + --- + ln)x[' ■■■X^- = J2k ■■■X^- = J2khkiXi,...,Xn), (3.12) 

0</iH \-l„<p k=0 ^ln=k k=0 

N{p, n) can also be expressed via the complete symmetric functions, 

N{p, n) = T.l=okhk{x,,... x) ^ ^^^^^^ 

iT'p \X\^ • • • ; Xfi) i-j 

In order to further simplify (|3.13| ), note that according to (|3.8|) 

p-i 

hp— l(^Xi, . . . , Xni 1) ^ ^ h]^{X\^ ■ ■ ■ ■> 2^n+l)- 

k=0 



Hence, setting = 1 and using again (|3.8|) one has 

p~l p—l k 

. . . , x„, 1, 1) = ^ hk{xi, . . . , x„, 1) = ^ ^ /ig, 

fc=0 k=Q q=0 

where here and below hg = hg{xi, . . . ,Xn)- Therefore 

p— 1 k p—l p—l p—l 

/ip_i(xi,...,a;„,l,l) = ^^/ig = ^^/ig = (3.14) 

k=0 q=0 q=0 k=q (J=0 

Prom here and ( |3.8| ) we deduce 

p p—l p 



php{xi, ...,Xn,l)- hp-i{xi, . . . 1, 1) = ^phk - - k)hk = ^khk{xi, . . . ,a;„). 

(3.15) 



A;=0 k=0 k=0 

Combining (|3.13|) with (|3.15| ) we finally obtain 



TTr/ \ ^p- 1(2^1; • • • ; Xn, 1, 1) 1 r\ 

N[p,n)=p — -— . (3.16) 

'''pl^'^lj • • • ) Xm i-j 

As expected, the average number of particles accommodated in the system cannot exceed p. 
Similarly for the average energy E[p,n) of the system one has 

E{p,n)= V {e,h + --- + ej y'l) f " , 

and therefore 

n n 

E{p,n) = "^eiXid^^ \nhp{xi, . . . ,x„, 1) = ^e^Xj^^^ In Z{p,n). (3.17) 

i=l 1=1 

Let us determine the equilibrium distribution of the particles on an arbitrarily chosen 
orbital i. According to (|3.1(]| ), V{p, n;li, . . . , In) yields the probability for the system to be 
in the state (/i, . . . , /„), which means that /i particles are accommodated on the first orbital, 
I2 particles on the second, and so on. Therefore the probability V{p,n;li) that k particles 
are accommodated on the i-th orbital is 

E (3.13) 



For the average number of particles /j on the i-th orbital we have 
k = ^li'P{p,n;li) = k — 



Zip, n) 

li=0 0<h+- + ln<p ' 

' 0<ii+-+i„<p 
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Hence 

k = Xid^^ ln/ip(xi, . . . ,x„, 1) = Xid^^ \nZ{p,n) = rd^^ In Z{p,n), i = l,...,n. (3.19) 
It follows that the average number of particles Na{p, n) on, say, the first s orbitals is 

s s s 

Na{p, n) = J2^i = Y, ^i^-. In n) = ^ rd^^ In n). (3.20) 

i=l i=l 1=1 

Evidently, the average energy Ei of the particles on the i-th orbital is : 

Ei = eiXid^A-n.Z{p,n) = eiTdf,AnZ{p,n), 2 = l,...,n. (3.21) 

Let us note that the expression for the probability ( p.l8| ) can be also written in a more 
compact form : 

V{p, n; li) = ^ ^ Xi ■ ■ ■ x^_-^ Xj^^-^ ■ ■ ■ 



1 -^i-l-^j+l -^n 



' k=0 h+---+li^i+li+i+---+l„=k 

X^' 

777 7 / ^ h/^i^Xi, . . . , . . . , Xn). 

Z[p,n^ ^ — ^ 



fc=0 



Applying ( pj.Sj ) to the rhs, we obtain the required expression for the probability to have k 
particles accommodated on the i-th orbital : 

^ / . 7 \ ^p—li ("^1 5 ■ ■ ■ 5 '^i— 1 5 -^i+l • • • ■^n)'^i hp—i^ {Xl ; . . • 1 Xi—\^ 1, Xi-\.\ . . . Xn}X^ 



Z{p,n) hp{xi,. . . ,Xn,l) 

(3.22) 

Some other thermodynamical functions can be determined too. For instance, from the 
general expression for the entropy 

S(p^ n) = -^(P'^)-^r=iM ^ ^^^^ ^3 23) 



and ( |3.17| ) there comes : 

^ n k ^ 

S{p,n) = — ^(ej - iJ,i)li + kB\nZ{n,p) = — ^(ci - fJ.i)xid^^ \nZ{p,n) + kBlnZ{n,p), 

i=l ^ i=l 

(3.24) 

which can also be written as 

S{p,n) = keirdr + l) lnZ{p,n) = kBd^rln Z{p,n), (3.25) 
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or equivalently 

S{p, n) = kerdr In n) - — fi, (3.26) 

T 

where 

fi = -rlnZ(p,n) (3.27) 

is the thermodynamical potential, another relevant thermodynamical function (in order to 
be consistent with the notation used so far we have replaced in (|3.24|) - (|3.27|) the Kelvin tem- 



perature T with the fundamental temperature r = /c^T, being the Boltzmann constant). 
Before proceeding further with some particular cases of the Hamiltonian ( |2.13|) we make 



a small deviation in order to draw a parallel between A-statistics and Bose statistics. To 
this end introduce new creation and annihilation operators 

B{p)f = ^, 2 = l,...,n, peN, (3.28) 

in W{j),n). It is easy to verify that for large values of p these operators satisfy "almost 
Bose" commutations relations lITBi : 



[B{p)j ,B{p)^] = [B{p)^,B{p)j] = 0, exact commutators, (3.29) 
[5(p)r, 5(p)+] ~ if + /2 + . . . + Z„ < p. (3.30) 

Therefore the representations of B{p)f in the Fock spaces W{p,n) with large values of p, 
restricted to states with a small number I1 + I2 + ■ ■ ■ + ln <^ p of accommodated particles pro- 
vide good approximations to Bose creation and annihilation operators (in finite-dimensional 
spaces). For this reason the operators B{p)f are said to be quasi-Bose creation and anni- 
hilation operators (of order p). In the limit p — >■ 00 these operators become indeed Bose 
operators [|16|. Therefore, parallel to quon statistics (see p4[ and references therein), A- 



statistics (for large values of p) can be considered as a theory allowing small violations of 
canonical quantum statistics in nonrelativistic quantum field theory. 

Coming back to the macroscopic considerations, we observe that for t = 1 the rhs of (|3.7|) 
reduces to the Bose GPF Zsosein) of a system with n orbitals, which are filled independently 
of each other : 



00 

fc=0 



Therefore, see (0| 



hk{xi, ...,Xn) = -— — . . . -— T = ZBosein). (3.31) 

yL XijyL X2) • • • yi Xn) 



ZBosein) - Z{p,n) = ^ hk{xi, . . . ,Xn). (3.32) 

k=p+l 

For sufficiently large values of p the rhs of (|3.32|) , which is always positive, can be made 
smaller than any positive number and therefore can be neglected. This is another confirma- 
tion (now from a macroscopic point of view) that A-statistics reduces to Bose statistics as 
the order of statistics p becomes large. An analogue of equation (|3.32|) for g-statistics is also 
available |]25|, eq. (5.14)]. 



In the next sections we shall consider some examples, the first one with all energies equal 
to each other. 
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4 The most degenerate case 

Here we consider an ensemble of particles with a Hamiltonian 



i7 = e^iV„ (4.1) 



1=1 



i.e. all orbitals have the same energy, and additionally we assume that they all have the same 
chemical potential, i.e., 

ei = €2 = . . . = e„ = e, /i^ = . . . = /i„ = /i =^ xi = X2 = . . . = Xn = x. (4.2) 

In this case the orbitals label internal degrees of freedom of the particles (spin, color, flavor) 
or, as more particular examples, the local orbitals of any multi-band Hubbard model or 
SU{N) Heisenberg chain. 

Most of the thermodynamical functions follow directly from the results of the previous 
section after the specialization (|4.2| ), but they can be written in a more explicit form. To 
this end one has to take into account that the number of terms in the rhs of ( |3.5| ) is {k + 
n — 1)!/A;!(ri — 1)!. Therefore 

k + n- 1\ ^ 




^ ,x\ (4.3) 



Then equation ( |3.9|) yields 



Z(p,n) = /ip(x^_^^^, 1) = ^ 1^^^^ 



x\ (4.4) 



k=0 

n times 



This sum can be rewritten as 



^(^'P) = E k r 5Z k 

k=0 ^ ^ k=p+l ^ ^ 



X 



1 fn+p\ fl,n + p + l > 

'x^^ 2Fi[ , ^ ; a; , (4.5) 



(1 - x)" \P+lJ V P + 2 



where 2-^1 is the classical hypergeometric function p6[. Compared to ([4.4| ), the expression 
in ( |4.5| ) looks at first sight a more complicated way of rewriting Z{n,p). Note, however, that 
the first term in the rhs of (14.51) is the Bose GPF 



^ 'k + n-l\ ,. 1 



fc=0 



^WBose = E( k y=0^- ^^-^^ 



Therefore, the second term is responsible for the difference between Bose and A-statistics. 
It carries, so to speak, the statistical interaction between the particles. 
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Using Euler's transformation formula for hypergeometric functions p6|, (1,3,15)], i.e. 

.F.f"''i.U(l-xr»-vJ^-"'^-';xV (4.7) 



c 

B~5|) can also be rewritten as 



The hypergeometric series appearing in ( [4.8| ) has the advantage that it is a terminating series 
(consisting of n terms), since one of its numerator parameters, 1 — n, is a negative integer. 
More explicitly, we can rewrite ( ^.8| ) as 



Equation ( [4.4| ) is convenient to deal with in those cases that the order of statistics p is a 
small number (and any number of orbitals n) . On the contrary, the expression ( [1.9[ ) is more 
appropriate for a relatively small number of orbitals (and any order of statistics p) . 
In the case of only one orbital, i.e. for the sl{2) GPF, equation ( [4.4|) yields 



p 1 p+i 

Z(p,l) = ^a;^' = — (4.10) 

k=0 



For the s/(3) GPF, the expression is 

2) = ±(k + 1).^ = + + P + (4.11) 



fc=0 



and it can be related to the Z{p, 1) partition function by 

This result can be further generalized. The GPF of sl{n + l) for any n can be related to the 
GPF of sl{2) : 

or equivalently 

^(^''') = (^T^ 2^ = (^T^ i-x ■ ^^-^^^ 

Clearly, after the specialization ( ^.2| ) the expression ( p.ll| ) for the average number of 
particles reads 

N{p,n) = xd^ln Z{n,p) (4.15) 
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and E{p,n) = eN{p,n). Another expression follows from ( |3.13|) , (f4.3| ) and ( [4.4| 



N{p,n) = ^'=' J . (4.16) 



Using the definitions of hypergeometric functions, (|4.16| ) can be rewritten as 



7^-^ -(p+1) 2F1 f ; X 

X 



Applying Euler's transformation to each of the 2-F1 functions yields an expression with ter- 
minating hypergeometric series in the numerator and denominator : 



N{p,n) = A" ^^x - (4.18) 

p+l,l-n. ^ 

p+2 ' 



So we find 



p!-(n+rt!ELo(-l)' 



'^(p.") - ;v v;tr ' ^^g^ i ■ (4.19) 

l-X - (^ + P)! Efc=o(-l)' (p+fc+l)fc!(n^fc-l)! , 

The last expression for N{p,n) is more appropriate to work with for small values of n, 
whereas ( [4.16| ) is more suitable for small values of p. 



From (|]2|) and we can also compute the probabiliy V{p, n; U) that U particles are 
accomodated on the i-th orbital : 

Then the average number of particles accommodated on the z-th orbital is 

As it should be, the result does not depend on the number i of the orbital : l\ = ■ ■ ■ = h = 

In — I- 

Using the binomial identity 
one verifies that the consistency condition N{p, n) = nl holds too. Hence 
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n. (4.23) 



Other thermodynamical functions follow straightforwardly. Equation ( p.24| ) for the entropy 
reduces to : 

S{p,n) = —{e~i2)N{p,n)-—n, (4.24) 
r r 

where Q = — t In Z{p,n) is the thermodjTiamical potential (|3.26 ). 

Let us consider in some more detail the dependence on the energy of the average number 

of particles in the system N{p,n), i.e. the distribution function. As an energy variable we 

take 

y = x = e-^ (4.25) 

namely the energy in units of r. We will consider two extreme cases. 
• n = 1 and any p : 

Note that at p = 1 one obtains the Fermi-Dirac distribution, 

N{IA) = —^, . (4.27) 

e^r + 1 



In Figure 1 we plot N{p, 1) for p=l,2, 3, 4, 6, 8, 12, 16, oo. 




The lowest curve {p = 1) yields the Fermi-Dirac distribution. Increasing p from 1 to 
oo one "deforms" it into the Bose-Einstein distribution [p = oo) 



N{p, oo) = 




(4.28) 
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• p = 1 and any n : 

— Tl 

N{l,n) = -^ . (4.29) 

N{l,n) is always smaller then 1, so the system can accommodate at most one par- 
ticle. As an example we plot the distribution functions for a system with n = 
1, 2, 4, 8, 16, 32, 64, 128 orbitals (Figure 2). 




The first curve n = 1 (from the left) corresponds to the Fermi-Dirac distribution 
function. With the increase of the number of orbitals the average occupation number 
of the system increases for fixed y. In particular for e = /i we have that iV(l, n) = 
All curves are "Fermi-like" but the half-filling is shifted to the right, at y = Inn. 

It should be noted that the curves in Figure 2 give the average number of particles in 
the system, not on a certain orbital. 

The particles described above [p = 1, n > 1) are called hard-core bosons. They appear 
naturally in various models of condensed matter physics and nuclear physics (for more 
discussions and references see |]TB[)- 



5 Equidistant energy levels 



Let us now consider the Hamiltonian p.l3| ) with equidistant energies e^. Denote the gap 



between the different energy levels by A > 0. This means that €2 = ei + A, €3 = ei + 2A, 



etc., or 



ei + («-l)A, (« = l,2,...,n). (5.1) 



Just as in the previous section, we shall assume that yUi = /i2 = ■ ■ ■ = /Xn = In this setting 
the different orbitals correspond to different energy levels. Following the notation of 
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we have 



Xi = exp 



11 -ei 



T 



exp 



/i - ei 



T 



A 

exp ( 



i-l 



xq 



J-i 



where we have used the notation 



X = Xi = exp 



/i - ei 



r 



and q = exp 



A 

r 



(5.2) 



(5.3) 



In order to write down the grand partition function, we can use ( p.9|) and the speciahza- 
tion given above : 

p 

Z{p, n) = hk{x, qx, q'^x, . . . , q"'~^x) = hp{x, qx, q'^x, . . . , q^~^x, 1). (5.4) 

fc=0 

The symmetric functions simphfy under this speciahzation. To see this, consider their gen- 
erating function ( p.7| ). Since [^, p. 26] 



1 



n + k — 1 
k 



xh\ 



{I - xt){l - qxt) ■ ■ ■ {I - q^-^xt) 

where ["^j denotes the g-binomial coefficient or Gaussian polynomial |2^, p. 26] : 

(1 - g™)(l - g™"i) ■ • ■ (1 - g™-*--+i) 



(5.5) 



m 
k 



(l-g)(l-g2)...(i_g/c) 



(5.6) 



it follows from (13.71) that 



hk{x,qx,q x, . . . , x) 



n + k — 1 
k 



x^. 



(5.7) 



Observe that in the limit g — 1, the g-binomial [™] goes to the ordinary binomial coefficient 
{X). Using (13), dp) implies 



Z(p,n) = ^ 



A:=0 



n + — 1 



(5. 



Using the g-raising factorials P7[ , 

(a; g)fc = (1 - a)(l - ga) ■ • ■ (1 - g^-^a), (5.9) 
and the classical g-generalized hypergeometric series, called basic generalized hypergeometric 



series 



27|, |26|] , this can be rewritten as 



=2^7 ^a; =2$i 



fc=0 



g", g P 



(5.10) 
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The average number of particles in the system follows from ( p.l3| ) : 

This expression cannot be further simplified. 

Another quantity that carries relevant information about the system is the average num- 
ber of particles accommodated on a particular orbital. Let li be this average for the i-th 
orbital, i = 1,2, ... ,n. Following ( |3.19| ), we have 

k = —^xAXZ{p,n)), (5.12) 
z [p, n) 

in which we have to substitute Xi = q^'^x. This expression can be written in the following 
more explicit form : 



n + I — I 



I 



xK (5.13) 



The derivation of ( ^.13|) , which is not so trivial, is given in the Appendix. 

The main conclusion from ( |5.13| ) is that the "population" of the orbitals depends essen- 
tially on their level i via where q = exp (—7-) < 1 : as i grows, the average number 
of particles U decreases. Otherwise said : the higher the energy level, the lower the average 
number of particles. 

If we consider the extreme case p = 1, where the system contains only one particle, and 
any n (the other extreme case, any p and n = 1 coincides with the most degenerate case) 
then there comes 

^ ' ' e/5(^i-^) + (l + g + --- + g"-i) ^ r ^ ' 

The case g = 1 (A = 0) corresponds to the degenerate case. 

For values of g = exp(— A/r) ^ 1, i.e., for large gaps between the energy levels or very 
low temperature, one can neglect all positive powers of q in ( p.l4| ). What remains is the 
Fermi-Dirac distribution 

Nil,n)^ ^, . (5.15) 

Continuing with this extreme case (where p = 1), the expression for the average number 
of particles on orbital i reads 

li = —T, ^ ; TT: i = \,...,n. (5.16) 

For very low temperatures, or big A, ( |5.16| ) reduces to 

Ix ^ -57 — \ and ^ if i> 1. (5.17) 

The latter means that if the system contains a particle, it is "sitting" permanently on the 
first, i.e. on the lowest energy orbital. This also explains why iV(l,n) ~ l\. 

The expressions for the entropy Sij), n) and the thermodynamical potential f2(p, n) follow 
from ( |3.25| )-( p.27| ) and cannot be simplified considerably. 
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6 Two species of particles 

We assume in this section that the system under consideration consists of two species of 
particles. Those of the first kind A (resp. of the second kind B) have one and the same 
energy and chemical potential fiA (resp. €3 and /i^)- To be more precise, the Hamiltonian 
of the system is 

m n 
1=1 i=m+l 

For convenience we consider a system with an even number of orbitals : n = 2m, m G N. 
The first m orbitals refer to single particle states of kind A, and the remaining m to single 
particle states of kind B. 

The probability for the system to be in a state r = (Zi, . . . , Z„) is given by (|0|) , which in 
this case reads : 

ll + ... + lmlm + l + --- + ln 

^^^■"■•^^ - ^-ZlKn) ■ 

where 

xa = exp ^ j , xb = exp ^ j . (6.3) 



In order to write down the grand partition function ( p.6|) we use the following identity 

k 

^fc('Cl) • • • ) •^rri) • • • ) •^n) ^ ^ hy{x\^ • • • ; -^m) ^fc— r (-^m+l ; • • • ; -^n); (^■'^) 

r=0 

which can easily be derived from the generating function ( |3.7| ). Then, in view of ( |4.3| ) 

fr + m — \\ fk — r + m — l\ 
hk[ XA, ■ .^.,Xa^ , Xb, . ,Xb) = } ^ ^ J. 



1" ^k- 



\x,x 



1-^2 



m times, m times ^ ° 



The latter can also be expressed by means of a hypergeometric function : 

fk + m — l\ ( rrij—k Xi . 
hkj xA, Xa , xb,..j, Xb) = \ ^ J ( 1 _ ^ _ ' ~ ) ^2- (6-5) 

m times, m times 
Hence the GPF (|3.6|) reduces to the following expression 

p k 

^ i r -\- 

r J \ k — r 



k=0 r=0 

or 



p 

xi (6.7) 



7f \ ^fk + m-l\ f m,-k xA 
= k )'^'[l-m-k^Vj 
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An immediate consequence of ( P . 1 1| ) is the expression for the average number of particles 
in the system : 

N{p,n) = {xaOxa + XBdxB)\nZ{p]n) = r^df,^ + df,g) In Z{p,n). (6.8) 
Using ( 3.13| ), (|6.5| ) and (6/7), one can write N{p,n) in a more exphcit form : 

Nip, n) = ) ^ . (6.9) 

Ep fk+m—l\ p / m.,—k . xb_ \ 
k=0{ k ^'2-^1 \^l-s-fc' xb) 



k 



From ( |6.2|) one derives the probabihty V{p,n; Ma, Mb) for the system to contain Ma 
particles of kind A and Mb particles of kind B : 

«; M.) ^ ^(-^ ;™ - ^) - ^).-.-. (e,.. 

Consequently 

yields the probability for the system to accommodate Ma particles of kind A. Therefore, 
the thermal average of the particles of kind A reads : 

Ma=0 Mb=0 ^ A / \ a / 

Formulas ( |6.11| ) and ( |6.12| ) can be re-expressed in terms of a hypergeometric function : 
V{p,n-MA) = r ^ (6.13) 



Z{p,n)\ Ma 

1 p-Ma+1 + P - Ma\ fl,m + p- Ma + 1_ \ \ Ma 



j=^'( p-M, + 2 

and 

More formally, we can also write 

N{p,n)A = XA<9^^1nZ(p;n) = Tdf,^\nZ{p,n). (6.15) 
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The termal averages E(p, n)A and E{p, 11)3 of the particles of kind A and B are evident 



and therefore 



E{p, n)A = tAN{p, n)A, E{p, 11)3 = esN^p, 71)3, 



E{p, n) = E{p, n)A + E{p, 71)3 



(6.16) 
(6.17) 



yields the average energy of the system. 

From ( p.24[ )-( p.27D the expression for the entropy follows : 



S{p,n) 



— {(^A - IJ-a)N{p, n)A + — (es - fiB)N{p, n)^ - — fi, 

T T T 



(6.18) 



where fl = —T\nZ{p,n) is the thermodynamical potential. 

It is instructive to consider an example in more detail. Let us fix n = 4 and also take 
p = 4. We shall draw a graph of the average number of particles N{p,n) = iV(4,4) (see 
equation (|6.9| )), as a function of two energy variables yA and associated with the two 
kinds of particles of the system, i.e. 



VA 



Vb 



(6.19) 



r r 

This graph is given in Figure 3. Clearly, this graph is symmetric with respect to yA and ys- 



Figure 3: Graph of A^(4, 4) for yA and j/s in the range [—5, 5]. 




Let us now also consider, for this same example, the graph of the average number of 
particles of kind A, i.e. iV^(4, 4). The expression follows from (|6.12|) . The graph is given in 
Figure 4. 

Comparing Figure 3 with Figure 4, one can make a distinction between four different 
regions in terms of the energy variables yA and ys- The sector {yA < 0,yB > 0) is populated 
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Figure 4: Graph of A^^(4, 4) for ua and ys in the range [—5,5]. 




mostly with particles of kind A, and the sector (y^ > 0,^^ < 0) mostly with particles of 
kind B. In the sector {yA < 0,yB < 0), the population of particles of kind A and of kind B 
is approximately the same. Finally, the sector {yA > 0,yB > 0) is essentially unpopulated. 
The average number of accommodated particles is never bigger that 4, as it should be, since 
p = 4. 



7 Concluding remarks 

In the present paper we have studied the thermal properties of "free" particles, which interact 
only via statistical interaction. The latter stems from the restrictions imposed by the Pauli 
principle : the system under consideration cannot accommodate more that p particles if the 
order of statistics is p. This property holds independently of the number of orbitals; there 
can even be infinitely many. 

By definition A-statistics is closely related to certain (more precisely, symmetric or Fock) 
representations of the Lie algebra sl{n + 1), including n = oo. Apart from that, A-statistics 



belongs to the class of exclusion statistics as defined in [^, Section 5]. Okubo |^ has also 



reformulated this in the language of Lie-triple systems. In |jT6| we have argued that under 



certain natural assumptions A-statistics can be interpreted as an exclusion statistics in the 
sense of Wu . 

Apart from the general case we have considered some specific examples. In particular we 
have shown that for n = 1 and any p the FD distribution function (n = 1, p = 1) deforms 
into the BE distribution function {n = 1, p = oo) with the growth of p, see Figure 1. In 
this case A-statistics reduces to Gentile statistics [|1| (see also |3^). In the more general case 
of any number of orbitals n the above picture is modified. In the limit p ^ oo one obtains 
again the Bose distribution function N{p = oo,n) = nx/{l — x). However at p = 1 the 
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distribution function N{p = l,n) is a distribution function of hard-core fermions, see ( [4.29|) , 
and not of fermions. 

Another observation to mention is in the case with equidistant energy levels. Without 
any input from quantum groups it turns out that the GPF is a g-deformation of the GPF 
of the most degenerate case. More precisely, the equidistant GPF ( p. 81 ) is obtained from the 
"nondeformed" GPF (|4.4| ) by a g-deformation of the binomial coefficients. Another property 
natural to expect, demonstrated here for p = 1, is that at very low temperatures the average 
number of particles of the system is the same as the average number of particles on the lowest 
energy level, which means that all allowed particles (in the general case p) "condensate" on 
the lowest level. 

Despite of the fact that A-statistics does not belong to the class of deformed Bose statis- 
tics, it yields a good approximation to Bose statistics. Apart from that the Fock spaces 
do not contain states with negative norm. Therefore, parallel to quons, A-statistics with 
large values of p is a good candidate for the description of small violations of Bose statis- 
tics in quantum field theory. Similarly as for quons however, we do not know how to 
satisfy the locality condition in relativistic quantum field theory. Therefore, one cannot ex- 



pect to derive relations between charge conjugation, unitarity and statistics as in [31]. It 



would be interesting to see whether such relations can be derived in the frame of causal 
A-statistics . 



Finally we point out that our considerations are incomplete in the sense of traditional 
thermodynamics, because we have not introduced the concept of volume and hence of pres- 
sure, etc. In our picture the volume can be introduced in several ways. One natural way 
would be to relate the order of statistics p to a unit volume V : if p is the maximal number 
of particles to be accommodated in V, then it is natural to assume that twice more particles 
could be accommodated in the volume 2V. This is one, but not the only plausible possibility. 
We shall return to this issue elsewhere. 



Appendix : Proof of equation ( |5.13| ) 

First, we wish to find an expression for hk{xi, . . . , Xj, . . . , a;„). The notation means that 
Xi has been removed from the list of variables ( x^, . . . , -£^n); SO hi^(^x\^ • • • 7 Xi^ • ' ' •} Xn] stands 
for a symmetric function in n — 1 variables. Multiplying ( p.7|) by (1 — Xjt), it follows easily 
that 

hk{Xi, . . . ,Xi,. . . ,Xn) = hk{x) - Xi/lfc_i(x), 

where hk(x.) = hk{xi,X2, ... ,Xn). 

Consider now the general expression for /j, as given in ( p.l9| ) : 

h = / Xid,,^Z{p,n), 
Z{p, n) 

or, using ( p.4|) , 

Xid^^Z{p,n) = 

0<h+l2 + - + ln<P 
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■^1 -^i-l-^i+l 



■ X. 



rf , , , ilf nf 

X1 ^i-l^iJ^l 



ii=0 0<h+-+k-i+k+i+-+l„<p 
P p-h 

= E'--?E E 

li=0 k=0 ih+---+li-i+k+i+-+l„=p-k) 

P p-h 

^ ^ ^i^^j ^ ^ hk{^X\^ ■ ■ ■ 1 Xi, • • • ; 3^n) 

ij=0 fc=0 

= ^ /jo;^' ^(/iA,.(x) - Xihk-iix)). 

In this last expression, we can make the speciahzation Xi = q^~^x. From ( ^.71) we know 
aheady how the functions hk{^) speciahze, so there comes (replacing also the summation 
variable U by /) 

p p-i 



E'('j"'-)'E 



1=0 



k=0 



n + k — 1 
k 



x^ - ^x 



n + k-2 
k-1 



Replacing g* by a new variable a, this can be rewritten as 



p p-i 





n + k — 1 




'n + k- 2' 


( 




— a 






k 





EE'"' 

1=0 k=0 

Collecting equal powers of a, this reduces to 

P P-l r 

E"'E 

1=1 k=0 

Putting back a = gives the relation (3.69), which we wanted to prove. 
Observe that one summation can be performed in (|A.1| ) : 



n + k — 1 
k 



X 



k+l 



(A.l) 



p p-i 

E"'E 

1=1 k=0 



n + k — 1 
k 



X 



k+l 



p-l 

E 

fc=0 

E 

k=0 



n + k — 1 
k 

n + k — 1 
k 



p—k 

E( 

1=1 



x^ > (ax)' 



X 



^,ax — [ax) 



\ — ax 



-)• 



Replacing again a by ^ yields an alternative expression for (|5.13|) . 
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